Predicting the Tolerated Sequences for Proteins and Protein Interfaces Using RosettaBackrub Flexible Backbone Design

Predicting the set of sequences that are tolerated by a protein or protein interface, while maintaining a desired function, is useful for characterizing protein interaction specificity and for computationally designing sequence libraries to engineer proteins with new functions. Here we provide a general method, a detailed set of protocols, and several benchmarks and analyses for estimating tolerated sequences using flexible backbone protein design implemented in the Rosetta molecular modeling software suite. The input to the method is at least one experimentally determined three-dimensional protein structure or high-quality model. The starting structure(s) are expanded or refined into a conformational ensemble using Monte Carlo simulations consisting of backrub backbone and side chain moves in Rosetta. The method then uses a combination of simulated annealing and genetic algorithm optimization methods to enrich for low-energy sequences for the individual members of the ensemble. To emphasize certain functional requirements (e.g. forming a binding interface), interactions between and within parts of the structure (e.g. domains) can be reweighted in the scoring function. Results from each backbone structure are merged together to create a single estimate for the tolerated sequence space. We provide an extensive description of the protocol and its parameters, all source code, example analysis scripts and three tests applying this method to finding sequences predicted to stabilize proteins or protein interfaces. The generality of this method makes many other applications possible, for example stabilizing interactions with small molecules, DNA, or RNA. Through the use of within-domain reweighting and/or multistate design, it may also be possible to use this method to find sequences that stabilize particular protein conformations or binding interactions over others

Quantifying the Ocean, Freshwater and Human Effects on Year-to-Year Variability of One-Sea-Winter Atlantic Salmon Angled in Multiple Norwegian Rivers

Many Atlantic salmon, Salmo salar, populations are decreasing throughout the species' distributional range probably due to several factors acting in concert. A number of studies have documented the influence of freshwater and ocean conditions, climate variability and human impacts resulting from impoundment and aquaculture. However, most previous research has focused on analyzing single or only a few populations, and quantified isolated effects rather than handling multiple factors in conjunction. By using a multi-river mixed-effects model we estimated the effects of oceanic and river conditions, as well as human impacts, on year-to-year and between-river variability across 60 time series of recreational catch of one-sea-winter salmon (grilse) from Norwegian rivers over 29 years (1979–2007). Warm coastal temperatures at the time of smolt entrance into the sea and increased water discharge during upstream migration of mature fish were associated with higher rod catches of grilse. When hydropower stations were present in the course of the river systems the strength of the relationship with runoff was reduced. Catches of grilse in the river increased significantly following the reduction of the harvesting of this life-stage at sea. However, an average decreasing temporal trend was still detected and appeared to be stronger in the presence of salmon farms on the migration route of smolts in coastal/fjord areas. These results suggest that both ocean and freshwater conditions in conjunction with various human impacts contribute to shape interannual fluctuations and between-river variability of wild Atlantic salmon in Norwegian rivers. Current global change altering coastal temperature and water flow patterns might have implications for future grilse catches, moreover, positioning of aquaculture facilities as well as the implementation of hydropower schemes or other encroachments should be made with care when implementing management actions and searching for solutions to conserve this species

Accessing ns–μs side chain dynamics in ubiquitin with methyl RDCs

This study presents the first application of the model-free analysis (MFA) (Meiler in J Am Chem Soc 123:6098–6107, 2001; Lakomek in J Biomol NMR 34:101–115, 2006) to methyl group RDCs measured in 13 different alignment media in order to describe their supra-τc dynamics in ubiquitin. Our results indicate that methyl groups vary from rigid to very mobile with good correlation to residue type, distance to backbone and solvent exposure, and that considerable additional dynamics are effective at rates slower than the correlation time τc. In fact, the average amplitude of motion expressed in terms of order parameters S2 associated with the supra-τc window brings evidence to the existence of fluctuations contributing as much additional mobility as those already present in the faster ps-ns time scale measured from relaxation data. Comparison to previous results on ubiquitin demonstrates that the RDC-derived order parameters are dominated both by rotameric interconversions and faster libration-type motions around equilibrium positions. They match best with those derived from a combined J-coupling and residual dipolar coupling approach (Chou in J Am Chem Soc 125:8959–8966, 2003) taking backbone motion into account. In order to appreciate the dynamic scale of side chains over the entire protein, the methyl group order parameters are compared to existing dynamic ensembles of ubiquitin. Of those recently published, the broadest one, namely the EROS ensemble (Lange in Science 320:1471–1475, 2008), fits the collection of methyl group order parameters presented here best. Last, we used the MFA-derived averaged spherical harmonics to perform highly-parameterized rotameric searches of the side chains conformation and find expanded rotamer distributions with excellent fit to our data. These rotamer distributions suggest the presence of concerted motions along the side chains

What explains the North–South divide in Italian tax compliance? An experimental analysis

This is the author accepted manuscript. The final version is available from Palgrave Macmillan via the DOI in this recordI undertake a comparative study assessing the North–South divide in Italian tax compliance, employing the largest behavioral tax compliance experiment to date. Contrary to a large body of literature, I argue that willingness to pay taxes is constructed within a specific institutional environment and reflects the country’s quality of institutions. To test this hypothesis, I use controlled tax compliance experiments from four laboratories in Capua, Rome, Bologna, and Milan. By employing the experimental method, I am able to hold institutions constant allowing me to isolate cultural variation. Contrary to cultural explanations for tax compliance, when controlling the institutional environment, there is no difference in tax compliance. Furthermore, using social value orientation to compare prosociality, I also find no differences between the two regions. I therefore conclude that individuals’ relationship to their states shapes these behavioral differences in tax compliance.Funds for this research were provided by the European Research Council (Grant Agreement No. 295675 )

In Silico Elucidation of the Recognition Dynamics of Ubiquitin

Elucidation of the mechanism of biomacromolecular recognition events has been a topic of intense interest over the past century. The inherent dynamic nature of both protein and ligand molecules along with the continuous reshaping of the energy landscape during the binding process renders it difficult to characterize this process at atomic detail. Here, we investigate the recognition dynamics of ubiquitin via microsecond all-atom molecular dynamics simulation providing both thermodynamic and kinetic information. The high-level of consistency found with respect to experimental NMR data lends support to the accuracy of the in silico representation of the conformational substates and their interconversions of free ubiquitin. Using an energy-based reweighting approach, the statistical distribution of conformational states of ubiquitin is monitored as a function of the distance between ubiquitin and its binding partner Hrs-UIM. It is found that extensive and dense sampling of conformational space afforded by the µs MD trajectory is essential for the elucidation of the binding mechanism as is Boltzmann sampling, overcoming inherent limitations of sparsely sampled empirical ensembles. The results reveal a population redistribution mechanism that takes effect when the ligand is at intermediate range of 1–2 nm from ubiquitin. This mechanism, which may be depicted as a superposition of the conformational selection and induced fit mechanisms, also applies to other binding partners of ubiquitin, such as the GGA3 GAT domain

Probing the Flexibility of Large Conformational Changes in Protein Structures through Local Perturbations

Protein conformational changes and dynamic behavior are fundamental for such processes as catalysis, regulation, and substrate recognition. Although protein dynamics have been successfully explored in computer simulation, there is an intermediate-scale of motions that has proven difficult to simulate—the motion of individual segments or domains that move independently of the body the protein. Here, we introduce a molecular-dynamics perturbation method, the Rotamerically Induced Perturbation (RIP), which can generate large, coherent motions of structural elements in picoseconds by applying large torsional perturbations to individual sidechains. Despite the large-scale motions, secondary structure elements remain intact without the need for applying backbone positional restraints. Owing to its computational efficiency, RIP can be applied to every residue in a protein, producing a global map of deformability. This map is remarkably sparse, with the dominant sites of deformation generally found on the protein surface. The global map can be used to identify loops and helices that are less tightly bound to the protein and thus are likely sites of dynamic modulation that may have important functional consequences. Additionally, they identify individual residues that have the potential to drive large-scale coherent conformational change. Applying RIP to two well-studied proteins, Dihdydrofolate Reductase and Triosephosphate Isomerase, which possess functionally-relevant mobile loops that fluctuate on the microsecond/millisecond timescale, the RIP deformation map identifies and recapitulates the flexibility of these elements. In contrast, the RIP deformation map of α-lytic protease, a kinetically stable protein, results in a map with no significant deformations. In the N-terminal domain of HSP90, the RIP deformation map clearly identifies the ligand-binding lid as a highly flexible region capable of large conformational changes. In the Estrogen Receptor ligand-binding domain, the RIP deformation map is quite sparse except for one large conformational change involving Helix-12, which is the structural element that allosterically links ligand binding to receptor activation. RIP analysis has the potential to discover sites of functional conformational changes and the linchpin residues critical in determining these conformational states

Using Entropy Maximization to Understand the Determinants of Structural Dynamics beyond Native Contact Topology

Comparison of elastic network model predictions with experimental data has provided important insights on the dominant role of the network of inter-residue contacts in defining the global dynamics of proteins. Most of these studies have focused on interpreting the mean-square fluctuations of residues, or deriving the most collective, or softest, modes of motions that are known to be insensitive to structural and energetic details. However, with increasing structural data, we are in a position to perform a more critical assessment of the structure-dynamics relations in proteins, and gain a deeper understanding of the major determinants of not only the mean-square fluctuations and lowest frequency modes, but the covariance or the cross-correlations between residue fluctuations and the shapes of higher modes. A systematic study of a large set of NMR-determined proteins is analyzed using a novel method based on entropy maximization to demonstrate that the next level of refinement in the elastic network model description of proteins ought to take into consideration properties such as contact order (or sequential separation between contacting residues) and the secondary structure types of the interacting residues, whereas the types of amino acids do not play a critical role. Most importantly, an optimal description of observed cross-correlations requires the inclusion of destabilizing, as opposed to exclusively stabilizing, interactions, stipulating the functional significance of local frustration in imparting native-like dynamics. This study provides us with a deeper understanding of the structural basis of experimentally observed behavior, and opens the way to the development of more accurate models for exploring protein dynamics

Tradeoff Between Stability and Multispecificity in the Design of Promiscuous Proteins

Natural proteins often partake in several highly specific protein-protein interactions. They are thus subject to multiple opposing forces during evolutionary selection. To be functional, such multispecific proteins need to be stable in complex with each interaction partner, and, at the same time, to maintain affinity toward all partners. How is this multispecificity acquired through natural evolution? To answer this compelling question, we study a prototypical multispecific protein, calmodulin (CaM), which has evolved to interact with hundreds of target proteins. Starting from high-resolution structures of sixteen CaM-target complexes, we employ state-of-the-art computational methods to predict a hundred CaM sequences best suited for interaction with each individual CaM target. Then, we design CaM sequences most compatible with each possible combination of two, three, and all sixteen targets simultaneously, producing almost 70,000 low energy CaM sequences. By comparing these sequences and their energies, we gain insight into how nature has managed to find the compromise between the need for favorable interaction energies and the need for multispecificity. We observe that designing for more partners simultaneously yields CaM sequences that better match natural sequence profiles, thus emphasizing the importance of such strategies in nature. Furthermore, we show that the CaM binding interface can be nicely partitioned into positions that are critical for the affinity of all CaM-target complexes and those that are molded to provide interaction specificity. We reveal several basic categories of sequence-level tradeoffs that enable the compromise necessary for the promiscuity of this protein. We also thoroughly quantify the tradeoff between interaction energetics and multispecificity and find that facilitating seemingly competing interactions requires only a small deviation from optimal energies. We conclude that multispecific proteins have been subjected to a rigorous optimization process that has fine-tuned their sequences for interactions with a precise set of targets, thus conferring their multiple cellular functions

Dynamism in the solar core

Recent results of a mixed shell model heated asymmetrically by transient increases in nuclear burning indicate the transient generation of small hot spots inside the Sun somewhere between 0.1 and 0.2 solar radii. These hot bubbles are followed by a nonlinear differential equation system with finite amplitude non-homologous perturbations which is solved in a solar model. Our results show the possibility of a direct connection between the dynamic phenomena of the solar core and the atmospheric activity. Namely, an initial heating about DQ_0 ~ 10^{31}-10^{37} ergs can be enough for a bubble to reach the outer convective zone. Our calculations show that a hot bubble can arrive into subphotospheric regions with DQ_final ~ 10^{28} - 10^{34} ergs with a high speed, up to 10 km s-1, approaching the local sound speed. We point out that the developing sonic boom transforms the shock front into accelerated particle beam injected upwards into the top of loop carried out by the hot bubble above its forefront traveling from the solar interior. As a result, a new perspective arises to explain flare energetics. We show that the particle beams generated by energetic deep-origin hot bubbles in the subphotospheric layers have masses, energies, and chemical compositions in the observed range of solar chromospheric and coronal flares. It is shown how the emergence of a hot bubble into subphotospheric regions offers a natural mechanism that can generate both the eruption leading to the flare and the observed coronal magnetic topology for reconnection. We show a list of long-standing problems of solar physics that our model explains. We present some predictions for observations, some of which are planned to be realized in the near future.Comment: 44 pages, 20 figure

The Long-Baseline Neutrino Experiment: Exploring Fundamental Symmetries of the Universe

Major update of previous version. This is the reference document for LBNE science program and current status. Chapters 1, 3, and 9 provide a comprehensive overview of LBNE's scientific objectives, its place in the landscape of neutrino physics worldwide, the technologies it will incorporate and the capabilities it will possess. 288 pages, 116 figuresMajor update of previous version. This is the reference document for LBNE science program and current status. Chapters 1, 3, and 9 provide a comprehensive overview of LBNE's scientific objectives, its place in the landscape of neutrino physics worldwide, the technologies it will incorporate and the capabilities it will possess. 288 pages, 116 figuresThe preponderance of matter over antimatter in the early Universe, the dynamics of the supernova bursts that produced the heavy elements necessary for life and whether protons eventually decay --- these mysteries at the forefront of particle physics and astrophysics are key to understanding the early evolution of our Universe, its current state and its eventual fate. The Long-Baseline Neutrino Experiment (LBNE) represents an extensively developed plan for a world-class experiment dedicated to addressing these questions. LBNE is conceived around three central components: (1) a new, high-intensity neutrino source generated from a megawatt-class proton accelerator at Fermi National Accelerator Laboratory, (2) a near neutrino detector just downstream of the source, and (3) a massive liquid argon time-projection chamber deployed as a far detector deep underground at the Sanford Underground Research Facility. This facility, located at the site of the former Homestake Mine in Lead, South Dakota, is approximately 1,300 km from the neutrino source at Fermilab -- a distance (baseline) that delivers optimal sensitivity to neutrino charge-parity symmetry violation and mass ordering effects. This ambitious yet cost-effective design incorporates scalability and flexibility and can accommodate a variety of upgrades and contributions. With its exceptional combination of experimental configuration, technical capabilities, and potential for transformative discoveries, LBNE promises to be a vital facility for the field of particle physics worldwide, providing physicists from around the globe with opportunities to collaborate in a twenty to thirty year program of exciting science. In this document we provide a comprehensive overview of LBNE's scientific objectives, its place in the landscape of neutrino physics worldwide, the technologies it will incorporate and the capabilities it will possess